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Abstract 



We present a simple feedback description of parametric instabilities which can be applied to a variety of optical systems. 
Parametric instabilities are of particular interest to the field of gravitational- wave interfcrometry where high mechanical 
quality factors and a large amount of stored optical power have the potential for instability. In our use of Advanced 
LIGO as an example application, we find that parametric instabilities, if left unaddressed, present a potential threat to 
the stability of high-power operation. 
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O 1- Introduction 

Though unseen in the currently operating first gen- 
nI eration interferometric gravitational- wave antennae (e.g., 
1—1 GEO [1], TAMA [2], Virgo [3], LIGO [6]), designers of 
i— —i "second generation antennae may be faced with instabili- 
^ .ties which result from the transfer of optical energy stored 



in the detector's Fabry-Pcrot cavities to the mechanical 



modes of the mirrors which form these cavity. 
i 1 The field of gravitational wave interfcrometry was in- 
troduced to the concept of parametric instabilities (PI) 
couched in the language of quantum mechanics [4]. We 
present a classical framework for PI which uses the audio- 
sideband formalism to represent the optical response of 
.the interferometer [7, 14]. This in turn allows the for- 
O^l malism to be extended to multiple interferometer config- 
.urations without the need to rederive the relevant rela- 
t— I tionships; an activity that has consumed considerable re- 
sources [9, 10, 11, 17]. 

2. Parametric Instabilities 



> 



X 



The process which can lead to PI can be approached 
d 'as a classical feedback effect in which mechanical modes 

_ i 

of an optical system act on the electro-magnetic modes of 
the system via scattering, and the electro-magnetic modes 
act on the mechanical modes via radiation pressure (see 
figure 1). 

We start by considering a single mechanical mode of 
one optic in the interferometer, and computing the para- 
metric gain of that mode [12]. The resonant frequency of 
this mode determines the frequency of interest for the feed- 
back calculation. 1 The nascent excitation of this mechan- 
ical mode, possibly of thermal origin, begins the process 



1 In general, all of the modes of all of the optics should be consid- 
ered simultaneously, but with typical mode quality factors greater 
than 10 6 it is very unlikely that two mechanical modes will partici- 
pate significantly at the same frequency. 
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Figure 1: Parametric instabilities can be described as a classical 
feedback phenomenon. A given mechanical mode interacts with the 
pump field to scatter energy into higher order optical modes. This 
interaction is represented with £5} above. Its strength is given by the 
overlap integral of the mechanical mode, pump mode, and scattered 
field mode (B m ^ n in equation 2). While circulating in the interferom- 
eter, the scattered field interacts with the pump field and mechanical 
mode via radiation pressure, introducing the overlap integral Cg) a sec- 
ond time. 



by scattering light from the fundamental mode of the opti- 
cal system into higher order modes (HOM) . The resulting 
scattered field amplitudes are 



TP -—A W R 

^0 



(i) 



where Ao is the wavelength and E pump the amplitude, of 
the "pump field" at the optic's surface, A m is the am- 
plitude of the motion of the mechanical mode, and B m ^ n 
is its overlap coefficient with the n th optical mode. The 
overlap coefficient results from an overlap integral of basis 
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functions on the optic's surface 



fofn(u m ■ z) dr A 



(2) 



where u m is the displacement function of the mechanical 
mode, fo and f n are the field distribution functions for the 
pump field, typically a gaussian, and the n th HOM. Each 
of these basis functions is normalized such that 



] dr = V 



\fn\ 2 df ± = 1 



(3) 



(4) 



where V is the volume of the optic. 2 

The interferometer's response to the scattered field can 
be computed in the modal basis used to express the HOMs 
via the audio sideband formalism. The resulting trans- 
fer coefficients represent the gain and phase of the optical 
system to the scattered field as it travels from the optic's 
surface through the optical system and back. The optical 
mode amplitudes of the field which returns to the optic's 
surface is 



to the transfer function of the mechanical system at its 
resonance frequency u) m , 
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where M is the mass of the optic and Q m the quality factor 
of the mechanical resonance. The open-loop-gain of the PI 
feed back loop is therefore, 
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The parametric gain 1Z is the real part of the open- 
loop-gain 
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n=0 



with the usual implication of instability if lZ m > 1, and 
optical damping if lZ m < 0. 3 
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where G n is the transfer coefficient from a field leaving 
the optic's surface to a field incident on and then reflected 
from the same surface. The transfer coefficient G n is com- 
plex, representing the amplitude and phase of the opti- 
cal system's response at the mechanical mode frequency. 
Computation of G„ for an optical system is discussed in 
more detail in section 3. 

The scattered field which returns to the optic closes the 
PI feedback loop by generating radiation pressure on that 
surface with a spatial profile that has some overlap with 
the mechanical mode of interest. This radiation pressure 
force is given by 
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(6) 



where P = \E pum p\ 2 , and c is the speed of light. The factor 
of 2 appears since the force is generated by the fields as 
they reflect from the surface. 

Finally, the radiation pressure which couples into this 
mechanical mode adds to the source amplitude according 



3. Optical Transfer Coefficients 

When the pump field is phase modulated by the me- 
chanical mode of the optic, upper and lower scattering 
sidebands are produced. In general, these scattering side- 
bands will have different optical transfer coefficients in the 
interferometer; and G~. The combination of scatter- 
ing sidebands which leads to radiation pressure is 



f ft— (-<+* 



(10) 



(see appendix B for derivation). This section will describe 
a general method for computing G^. 

Given scattering matrices S„ which contain transfer 
coefficients for the n th HOM of the scattering sidebands 
from one point to the next in the optical system, we have 



Gt = ej (I - 



(11) 



where the basis vector e x is the x column of the identity 
matrix I, and e% is its transpose. The index x is used 
to select the field which reflects from optic of interest, as 
demonstrated in the following section. 4 

3.1. An Example: Fabry Perot Cavity 

As a simple and concrete example, we apply the above 
formalism to a Fabry Perot cavity (FPC) of length L. Fig- 
ure 2 shows the configuration and indices for each of the 
fields in the FPC. 



2 A more general normalization of u m which allows for non- 
uniform density would include the density function of the optic in 
the integral, with the result equal to the total mass. 



3 Appendix A relates this result to the results found in previous 
works. 

4 A general framework for constructing scattering matrices is de- 
scribed in [5] and will not be reproduced here. 
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Figure 2: A simple Fabry Perot cavity. In the example caculation, 
fields arc evaluated at each of the numbered circles. Field 4 is used 
to compute the optical gain of the scattered field produced by, and 
acting on, mirror B. 
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The scattering matrices for the upper and lower side- 
bands are 
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(12) 



where tA is the transmissivity for mirror A, r a and re are 
the reflectivities of the mirrors, and = e l (<K± w ™ L / c ) j s 
the propagation operator. The reflectivity and transmis- 
sivity used in §^ are amplitude values and are related to a 
mirror's power transmission by t = \[T and r = \/l — T. 
The propagation phase depends on the phase of the n th 
HOM, cf) n , and the scattering sideband frequency ±u> m . 

If we wish to evaluate lZ m for a mode of mirror B, we 
would use 



o 
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(13) 



to select its reflected field, number 4 in figure 2. To ar- 
rive at a numerical result for lZ m , we adopt the following 
parameters 



P = 1MW 
T A = 0.014 
L = 3994.5 m 



A = 1064 nm 
T B = 10~ 5 
M = 40 kg 



which are representative of an Advanced LIGO arm cavity 
operating at full power. For the moment, we will consider 
a single optical mode, the Hermitc-Gauss TEM11 mode, 
and a single mechanical mode 



Qm = 10 7 

B m ,HGii = 0.21 



oj m = 2tt x 29950 Hz 
Wi = -5.434 



both of which are shown in figure 3. 

For this set of values, we evaluate the parametric gain 



G hgu = 0.554 + ^2.72 
Grgix = -0.052 + i2.75 



Gff Gn = 0.502 + i0.03 



Figure 3: A mechanical mode of an Advanced LIGO test mass near 
30kHz and a HG TEM11 optical mode. For the mechanical mode, 
surface displacement amplitude normal to the surface, u m ■ z, is 
shown. For the optical mode, the basis function /jjgii amplitude is 
shown. In both cases, red is positive, blue is negative and green is 
zero. The X and Y-axes on both plots are in centimeters. 



to find it is small and negative. 

Allowing the mechanical mode frequency to artificially 
vary from 20 kHz to 50 kHz we can plot lZ m as a function 
of ui m . The result is shown in figure 4. The resonance at 
27.4 kHz is the upper scattering sideband of the TEM11 
mode which has negative parametric gain, indicating opti- 
cal damping. At 47.7 kHz the lower scattering sideband of 
the TEM11 mode resonates, this time resulting in positive 
feedback, but not enough to produce instability. 

Fabry Perot Cavity Gain Factors 
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Figure 4: Optical gains G^ G11 and G^j G11 , and parametric gain 
TZm,HGli are shown as a function of mechanical mode frequency. 
The resonance of the upper scattering sideband at 27.4 kHz has neg- 
ative parametric gain, indicating optical damping. The lower scat- 
tering sideband resonance at 47.7 kHz has positive gain, but does not 
produce instability. 

Thus far only one optical mode and only one mechan- 
ical mode have been considered. Extending the computa- 
tion to higher order optical modes requires that we elabo- 



3 



rate our expression for <j> n to include the Gouy phase. For 
an arbitrary Hcrmite-Gauss mode 



K = 0o - O n c 



(14) 



where <fio is the propagation phase of the TEMOO mode, 
and O n is the mode order of the n th HOM. 

Considering other mechanical modes is a matter of 
computing the mode shapes and frequencies for the mir- 
rors which make up the FPC; we use those of an Advanced 
LIGO test-mass. 5 The result of the full calculation of lZ m 
for all mechanical modes between 10 kHz and 90 kHz, in- 
cluding HOMs up to 9 th order, is shown in figure 5. 

Fabry Perot Cavity Parametric Gain 
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Figure 5: Parametric gains TZ m are shown for mechanical modes 
between 10 kHz and 90 kHz. Red circles mark modes with positive 
parametric gain, while green circles mark those with negative gain. 
This calculation uses HOMs up to 9 order, but does not include 
clipping losses, discussed in section 4. 



mode, contributions to 7Z m from modes of order O n > 4 
are limited, and modes with O n > 9 are insignificant. 

A more complete description of losses due to apertures 
includes diffraction effects, but this requires a more com- 
plex and interferometer dependent calculation. Even bet- 
ter is to use the eigenmodes of the full interferometer, and 
their associated losses, rather than the Hermite-Gauss ba- 
sis. This level of detail may not be rewarded, however, 
since modes which differ significantly from their Hermite- 
Gauss partners do so as a result of significant losses, which 
in turn make them irrelevant to PI. 

4-1- An Example: Advanced LIGO 




4. Clipping Losses 

Thus far we have ignored losses in the optical system. 
A lower limit on the losses is given by the loss of power due 
to the finite size of the optics, known as clipping loss. For 
low-order optical modes, these losses are usually insignif- 
icant by design, but losses can strongly impact the para- 
metric gain when the contribution from high-order optical 
modes is dominant. 

In optical systems such as gravitational-wave interfer- 
ometers, in which the beam size on the optics is made as 
large as possible without introducing significant loss in the 
TEMOO mode, high-order modes tend to fall off the cav- 
ity optics. Specifically, for an interferometer designed to 
have a few parts-per-million clipping losses for the TEMOO 



5 The Advanced LIGO test-mass mechanical modes used in this 
and the following examples are the result of finite clement modeling. 
For a discussion of numerical and analytic methods for calculating 
test-mass mechanical modes, see [16]. 



Figure 6: Fields in a power and signal recycled Fabry Perot Michel- 
son. This optical configuration is common to many of the 2 nd gen- 
eration gravitational-wave detectors. 

As a more interesting example, we apply the above 
formalism to an Advanced LIGO interferometer. Figure 6 
shows the layout of the optical system and the assignment 
of field evaluation points (FEPs). In this case care has 
been taken to minimize the number of FEPs and to follow 
each one with a propagation operation. In this way we can 
number the propagation distances L x according to their 
associated FEP, and the propagation operations become 

n± — p i(0n,x±w m L x /c) 

The scattering matrices for this interferometer can be 
split into a diagonal propagation matrix populated by Pn X , 
and a mirror matrix populated by reflectivity and trans- 
missivity coefficients (essentially one r and one t in each 
column), as follows 
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Losses can be included in the scattering matrix §^ 
most generally by allowing M„ to vary for each HOM. A 
simpler approach which is sufficient for clipping losses is to 
add a diagonal matrix C„ to equation 15 which effectively 
adds loss to each propagation step 



■"n n 
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T-n,2 



where 



tn . . 



(16) 



(17) 



is the amplitude transmission of the aperture associated 
with each propagation step. For this example, a 17 cm 
aperture is assumed for all of the optics, except the beam- 
splitter for which we assume a 13.3 cm aperture. 6 

In order to compute the parametric gain as a function 
of mechanical mode frequency, as in figure 4 for the FPC 
example, we use the same values as above 

P=1MW A = 1064nm 
M = 40kg Q m = l0 7 

and add the transmission of the new optics 

Tix = T IY = 0.014 T EX = T EY = 10" 5 
T PR = 0.03 T SR = 0.2 T BS = 0.5 



new lengths 



£{1,2,3,4} = 3994.5 m 
£{5,6} = 4.85 m £{7,8} = 4.9 m 
£{9,10} = 52.3m £{n,i2} = 50.6m 

and phases 

00,(1-8,11,12} = 00,(9,10} = 

0G, {1,2,3,4} = 2.72 0G,{5,6,7,8} = 

0G,{9,io} = 0.44 G ,{ii,i2} = 0.35 

which represent 156° of Gouy phase in the arms, 25° in 
the power recycling cavity and 20° in the signal recycling 
cavity. The results are similar to the FPC alone, with neg- 
ative gain near 27.4 kHz and positive gain near 47.7 kHz. 
The addition of the rest of the interferometer, however, 
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Figure 7: Comparison of 72. m hgii f° r a Fabry Perot cavity and 
Advanced LIGO are shown as a function of mechanical mode fre- 
quency. The Advanced LIGO computation is shown twice, the green 
curve includes clipping losses, but has perfectly matched arm cav- 
ities, the red curve adds a 0.1% mismatch between the arm cavity 
Gouy phases. The calculation is limited to the modes shown in fig- 
ure 3, with the mechanical mode frequency artificially adjusted to 
highlight the resonance near 47.7 kHz. 



leads to a narrow regions of high parametric gain for the 
HG11 mode, see figure 7. 

As noted in [17], realistically imperfect matching of in- 
terferometer optics can significantly change the parametric 
gains of the system. Our formalism can be made to repro- 
duce this result by allowing the Gouy phase in one arm of 
the interferometer to differ from that of the other arm, 

0n,{3,4} = (1 + e)0n,{l,2} 

where ey is the fractional departure of c/> n ,{3,4} from their 
nominal value. For our Advanced LIGO example, where 
radius of curvature errors of a few meters are expected, we 
take e = 10~ 3 . This tiny difference is sufficient to move 
the sharp features associated with a perfectly matched in- 
terferometer by more than their width, thereby changing 
the result of any PI calculation (see figure 7). These re- 
sults are similar to those in [8] , though their model of the 
power-recycling cavity was somewhat less general. 

Finally, the full calculation for Advanced LIGO is plot- 
ted in figure 8. To show the effect of clipping losses, modes 
which have lZ n > 1 in the absence of clipping are shown 
and connected to their clipped partners. 

5. Worst Case Analysis 

Evaluating the impact of PI on a gravitational- wave in- 
terferometer is complicated by the sensitivity of the result 
to small changes in the model parameters. In particular, 



6 The beam-splitter aperture includes the aperture presented by 
the electro-static actuators on IX and IY. 
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Figure 8: Parametric gain for all modes of an Advanced LIGO test- 
mass between 10 kHz and 90 kHz. To show the effect of clipping 
losses, a few of the modes which have significant gain in the absence 
of clipping are also included. Empty circles represent parametric 
gains computed without clipping. They are attached to their clipped 
partners, represented by filled circles, with dashed lines. 



uncertainty in the radii of optics used in the arm cavities 
lead to changes the Gouy phases which, while quite small, 
are sufficient to move the optical resonances of the cavity 
by more than their width. Similarly, mechanical mode fre- 
quencies produced analytically or by finite element model- 
ing may not match the real articles due to small variations 
in materials, assembly, and ambient temperature. This 
section describes a robust means of estimating the "worst 
case scenario" for a given interferometer. 

A simple approach to the "worst case" problem is to 
compute the parametric gain of each mode for multiple 
sets of plausible interferometer parameters. Varying all 
of the parameters is impractical and unnecessary as the 
results are primarily sensitive to the relative frequencies of 
the mechanical modes and the optical resonances in high- 
finesse cavities. 

In the case of Advanced LIGO, explored in section 4.1 
above, it is sufficient to vary the Gouy phases in the arm 
cavites by a 5 x 10 -3 , and the phases in the recycling 
cavities by a couple of degrees. We proceed with a Monte- 
Carlo type analysis in which we randomly vary the Gouy 
phases around their nominal values. We repeat the pro- 
cess for 120 thousand trials, then set an upper-limit on lZ n 
for each mode at the lowest value greater than 99% of the 
results (see figure 9). 7 This provides us with a trial num- 
ber insensitive statistic, the accuracy of which is limited 
primarily by the fidelity of our model. 

We find that Advanced LIGO faces the possibility of 
a few unstable modes. Among the 120 thousand cases 
considered, the mean value of the maximum 1Z„ among all 
modes was 5.8 and 99% of cases had a maximum paramet- 



ric gain value less than 45. Considering each mechanical 
mode independently, there are 32 modes which have the 
potential to be unstable and that all of the highly unsta- 
ble modes are between 15 kHz and 50 kHz. Taking into 
account that 4 test-masses make up the Advanced LIGO 
detector, we find the mean number of unstable modes to 
be 10 (2.5 per test-mass), with 99% of cases having 6 or 
fewer unstable modes per test-mass. 

One must keep in mind, however, that many of the 
parameters used in this model are adjustable (e.g., power 
level in the interferometer, mirror temperature and thus 
mechanical mode frequency and radius of curvature, etc.) 
and others are speculative (e.g., the quality factor of a me- 
chanical mode may depend strongly on its suspension [13, 
15]). Since the parametric gain scales directly with both 
power in the interferometer and with mechanical mode Q, 
we have chosen round values for these parameters which 
can be refined as higher fidelity numbers become known. 

Advanced LIGO, 99% Confidence Upper Limit 
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7 To speed the computation, only optical modes with some overlap, 
, and modest clipping losses, £? , > 0.7, arc considered. 



Figure 9: Worst case parametric gain for all modes of an Advanced 
LIGO test-mass between 10 kHz and 90 kHz. There are 32 potentially 
unstable modes, and more than 200 modes with TZ n > 0.1. 



6. Conclusions 

Parametric instabilities arc of particular interest to the 
field of gravitational-wave interfcrometry where high me- 
chanical quality factors and a large amount of stored opti- 
cal power has the potential for instability. We depart from 
previous work by constructing a flexible analysis frame- 
work which can be applied to a variety of optical systems. 
Though our examples use a Hermite-Gaussian modal basis 
to describe the optical fields, this formalism can be imple- 
mented using the modal basis best suited to the optical 
system at hand. 

In our use of Advanced LIGO as an example appli- 
cation, we find that parametric instabilities, if left unad- 
drcssed, present a potential threat to the stability of high- 
power operation. We hope that future work on solutions 
to parametric instabilities will be guided by these results. 



A. Comparison with Previous Works 

The notation used herein has been chosen to match the 
seminal work by Braginsky where possible. To relate this 
work to [4] . which considered a single Fabry-Perot cavity, 
wc need the following equalities 



2P = So 
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L(l + Auo 2 /S 2 ) 



B m,n - A 



M = m 



where symbols on the left side of each equation are defined 
herein, and the symbols on the right are used in their equa- 
tions 3, 4, 8 and 9. Braginsky uses to for the total mass 
of the optic, but here that symbol is used as the mechani- 
cal mode index, so we use M instead. For reference, their 
equation 4 is 

2£oQiQ m 



n 



mL 2 L 



such that our 1Z is the same as the left-hand side of their 
equation 8. That is, 



n A 



where the symbols on the right are those of Braginsky. 
One should note that Braginsky only considered one me- 
chanical mode and one optical mode, so his subscript 1 
replaces our to, n and the summation over modes is ab- 
sent. Similarly, B m „ in this work is the same as B in 
[12]. 

B. Detailed Derivation of F ra d 

We start by examining more carefully equation 1 in 
which the scattered field is generated by phase modulation 
of the pump field. A complete description of the pump field 
reflected from the optic is 



£pump{j '_!_) E pump fo{r J_ ) 
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(18) 



where fo(r±) is the normalized basis function of the pump 
field, and ujq = 2itc/\q is the frequency of the pump field. 
Mechanical motion of the reflecting surface, with a spatial 
profile and frequency defined by a resonant mode of the 
optic 

Z m (r±) = A m {u m {ri_) ■ z) e 1 ^ 1 (19) 

results in phase modulation of the pump field. The scat- 
tered field is expressed to first order as a sum of HOM 
fields 

oo 

e.cat(r±) = E scat , n f n (r ± ) e^e*"-* + e"^ 4 ) (20) 
with HOM basis functions f n (rx), and HOM amplitudes 

E scat,n ~~ 7 A m Ep Urn p]3 mn . (^1) 

Ao 



Note that in equation 20 there are 2 frequency compo- 
nents produced; these are the upper and lower scattering 
sidebands. 

The radiation pressure which couples the scattered field 
into mechanical motion is given by 

Prad(r±) = \ Srefi(r±)£; efl (r±) (22) 
where the reflected field is 

£refl(r±) = Spump(r±) + £rf rn ( r ±) + ^rtrn( r ±) ( 2 3) 

and are tnc return fields, with HOM amplitudes 

given by applying G„ to the scattered field. By defin- 
ing £^r trn as part of the reflected field, G^ represents the 
closed-loop gain of each scattering field. 

The component of radiation pressure which is at the 
mechanical mode frequency is 
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Integrating over the surface to get the force which couples 
into the mechanical mode of interest simplifies the notation 
somewhat to 

2 / 00 
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where the last step comes from substitution using equation 
21. We recover equation 6 by setting 
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